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£h . Abstract. Transport in disordered armchair graphene nanoribbons (AGR) with 

long-range correlation between quantum wire contact is investigated by transfer 
matrix combined with Landauer's formula. Metal-insulator transition is induced by 
disorder in neutral AGR. Thereinto, the conductance is one conductance quantum 
for metallic phase and exponentially decays otherwise, when the length of AGR 
. approaches to infinity and far longer than its width. Similar to the case of long- 

range disorder, the conductance of neutral AGR first increases and then decreases 
while the conductance of doped AGR monotonically decreases, as the disorder strength 
increases. In the presence of strong disorder, the conductivity depends monotonically 
and non-monotonically on the aspect ratio for heavily doped and slightly doped AGR 
respectively. For edge disordered graphene nanoribbon, the conductance increases with 
the disorder strength of long-range correlated disordered while no derealization exists, 
since the edge disorder induces localization. (Some figures in this article are in colour 
, only in electronic version) 
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1. Introduction 

Graphene, discovered in 2004 [JJ El El S], stimulated intensive research interest in 
fundamental properties and its potential application in nanoelectronic device due to 
its superior properties such as high mobility [3] and heat dissipation [5j. For an infinite 
pure graphene sheet, the electronic spectrum is linear around Dirac point, and can be 
described as massless Dirac fermions. Due to its unique electronic structure, many 
interesting phenomena had been observed such as half-integer quantum Hall effect 
[I] and bipolar current [6]. Now it is commonly thought that the most possible 
application may come from graphene nanoribbon, and the basic types of graphene 
nanoribbon are armchair-edged and zigzag-edged ones. Defects and impurities are 
inevitable in graphene-based materials, therefore effects of defects, impurities and 
disorder on transport properties are of application interest and had been reviewed in 
Refs. [3 El El HO]. As for disorder, it is usually classified as short- and long-range 
disorder depending on the ratio of the interaction range to the lattice distance [IT] . 
Anderson disorder is a typical short-range disorder, and localization is induced by 
strong disorder. Long-range disorder such as Guassian-type disorder [12J induces metal- 
insulator transition at Dirac point. After several years of intensive investigation, it 
was well known that the scattering between two nonequivalent valleys is absent in the 
presence of long-range disorder [9] . Interestingly, off-diagonal disorder [131 E] enhances 
electronic transport and superconductivity correlation at Dirac point. Recently, more 
and more attentions were paid to the effect of charged impurities on the transport 
in graphene based materials [151 US] since the observation of electron-hole puddles in 
graphene [17] . 

It was commonly believed that Anderson localization [18] in one- and two- 
dimensional systems is induced by even a weak disorder from the scaling theory [19] . It 
was found that one or several resonant states may occur in one-dimensional disordered 
chain with short-range correlation [201 Ell E2] . Delocalized states within a continuous 
energy range exist in one-dimensional system, when disorder in onsite energy is long- 
range correlated [231 El] , where the site energies have an approximate spectral density 
of the form s(k) oc k~ a with s(k) being the Fourier transformation of the two-point 
correlation function of site energies < e^ej > and k the inverse of the wavelength A of the 
undulation on the site-energy landscape. The scaling behavior of the resistivity verified 
that delocalized states persist in the thermodynamic limit, and the relation between 
disorder correlation and the resistivity proved that disorder correlation is important for 
derealization [25]. There exists a Kosterlitz-Thouless type metal- insulator transition 
in two dimensional systems with long-range correlated disorder [261 [27]. Nowadays, the 
effect of disorder with long-range correlation on transport is under intensive investigation 
in many fields, which recently had been reviewed in Ref. [28] . 

Recently, the effect of disorder with long-range correlation in graphene attracts 
many researchers' attention [29l [30l [311 32J. For example, transport in disordered 
graphene superlattice with long-range correlation was studied on basis of massless Dirac 
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equation [3T] . It was found that there exists a truly metallic phase at the thermodynamic 
limit in the presence of long-range correlation, and the incident angle of metallic state 
increases with the correlation strength. Furthermore the metallic phase was obtained 
in term of correlation strength and disorder strength. 

It was found that there exists equivalent mode selection between armchair graphene 
nanoribbon (AGR) and quantum wire (QW) [33]. However, there is no such equivalence 
between zigzag graphene nanoribbon (ZGR) and QW, and even-odd parity occurs when 
ZGR is connected to QW contact even though ZGR is metallic [3U [35]. Moreover, 
the conductance of graphene nanoribbon at Dirac point is independent of the contact 
structure, provided that there are dense modes around Dirac point [351 EE]. Similar 
to our previous work [351 EZ]? QW is adopted to mimic metallic contacts here. In 
this Letter, we investigate the transport of disordered AGR with long-range correlation 
based on tight binding model by transfer matrix method. Our main results are listed as 
follows. Localization-delocalization transition is induced by the long-range correlation 
of disorder when the disorder is weaker than one certain critical value. Thereinto, the 
conductance for metallic phase approaches one conductance quantum and exponentially 
decays otherwise, as the length increases to infinity and is far longer than the width. 
Similar to the case of long-range disorder [11] , the conductance of neutral AGR first 
increases and then decreases while the conductance of doped AGR monotonically 
decreases, as the disorder strength increases. For disordered AGR with moderate aspect 
ratio at strong disorder, /3 varies from negative to positive in the scaling function of 
conductivity 0(LP) outside and within a certain range of energy shift respectively. 
Therefore the conductivity depends monotonically and non-monotonically on the aspect 
ratio for heavily doped and slightly doped AGR respectively. 

The structure of this Letter is as follows. In Sec. II, we introduce the structure of 
AGR between QW contacts, the tight binding Hamiltonian and long-range correlated 
disorder; in Sec. Ill, we show the dependence of the conductance on the strength of 
correlation and disorder, and the dependence of the conductivity on the aspect ratio in 
the presence of strong disorder that induces Anderson localization at the thermodynamic 
limit; and in Sec. IV, we make a summary. 

2. Structure and formalism 

The geometry of AGR sandwiched by two normal metal contacts is illustrated in Fig. 1, 
where semi- infinite QW stands for the normal contact. The shape of AGR is determined 
by L and M carbon atoms in the x and y direction, respectively. Correspondingly, the 
length is L x = *Mi2l and the width is L y = (M — 0.5)a, where a = 2.46A is the lattice 
constant of graphene. Lattice sites in graphene are classified as two sublattices according 
to their topological structures with respect to their neighboring sites. 

The transport property of GNR is mainly contributed by tt orbital and the tight- 
binding approximation is valid to describe the electronic structure of GNR at the Fermi 
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Figure 1. Illustration of armchair graphene nanoribbon (AGR) between two semi- 
infinite quantum wire contacts. The red and blue circle stand for A and B sublattice 
in graphene ribbon respectively. Two inequivalent lattice vectors d[ and are also 
shown. The shape of GNR is determined by L and M carbon atoms in the x and y 
direction, respectively. Correspondingly, the length is L x = v ^ La and the width is 
L y = (M — 0.5)a with the lattice constant a being 2.46A. 



energy. The Hamilitonian is expressed by 

H = t ■ »u-/'K',A- (!) 

<ij,i'j'> ij 

where i and j are x- and y- index of lattice site respectively, the nearest neighbors hopping 
t is set as the energy unit, < ij,i'j' > denotes two nearest neighboring lattice site, /i 
is the chemical potential shift by the gate voltage applied to AGR and is zero at QW 
contacts, and Cjj (C^) are electron creation (annihilation) operator at the lattice site, 
respectively, is the onsite energy at lattice site (i, j), which is zero at QW contacts 
and disordered with long-range correlation at AGR. Two parameters a and 5 measure 
the strength of correlation and disorder respectively. As shown in Ref. [25], the onsite 
energy is defined as follows, 



N/2 

E 

k=l 



k' 



2tt 
iV 



l-a 



COS 



27rlijk 
N 



(2) 



Here N = L x M is the total lattice sites of AGR, and kj — (i — 1)M+ j is the sequencing 
of lattice sites that onsite energies at one column of lattice sites are as close as possible. 
The sequence of onsite energies is normalized with zero mean value < e >= and 
with the variance Ae = a/< e 2 > — < e > 2 = 5. The disordered onsite energies with 
long-range correlation are shown in Fig. |2j The disorder strength 5 is set as 1 and 
the parameter a determines the correlation strength between any two disordered onsite 
energies. The onsite energies are random as a is 1.0, and become smoother and finally 
similar to a sequence of stripped potential when a increases from 1.0 to 2.5. As shown 
in Fig. EJ^c), AGR is equivalent to a composite of several random p-n junctions. The 
hard-wall boundary condition is adopted outside of the system, and the conductance 
in multichannels' form is calculated as shown in Ref. [3H]. The Fermi energy E = is 
chosen as there are maximal channels at the band center of QW contacts. 
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Figure 2. On site energy in the presence of different long-range correlated disorder. 
The system size is 20 x 40. a is 1.0 (a), 2.0 (b) and 2.5 (c). The strength of disorder 
is 8 = 1. 



3. Numerical results and discussion 

3.1. The scaling behavior of the conductance at neutral point in the presence of 
long-range correlated disorder 
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Figure 3. Dependence of the conductance of neutral AGR on the length when a 
varies at S = 0.1 (a-b) and when 6 varies at a = 2.5 (c). (d) shows the conductance 
at /i=0 (square), 0.2 (circle) and 1.0 (upper triangle) with a = 2.5 and 5 — 0.1. The 
width of sample is 10. 



It is well known that the band gap of uniform AGR shows a three-family behavior 
from tight binding model and first principle calculation [39J, where these two methods 
predicted that AGR with N a = 3p + 2 as p being an integer has no band gap and finite 
band gap inversely proportional to the width respectively. Further it was proved that 
the band gap of AGR with N a = 3p + 2 approaches to zero when the width is wider than 
80 A through divide-and-conquer approach based on quasi-atomical minimal basis orbits 
(QUAMBOs) HQ]. iV a is 2M in AGR as shown in Fig. [[J and M = 3p + 1 stands for 
metallic AGR. Even for semiconducting AGRs with M ^ 3p + 1, the conductance and 
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the localization length always increases with a in the presence of long-range correlated 
disorder. For sake of simplicity, metallic uniform AGRs are chosen to investigate the 
effect of long-range correlated disorder here. 

Fig. [3] shows the dependence of the conductance on the length of AGR in 
the presence of several different correlated disorder. As shown in Fig. Efa-b), the 
conductance approaches to one conductance quantum at highly correlated weak disorder 
as the length increases, when a is larger than 2.0 at 5 = 0.1. This proves that long-range 
correlation in disordered onsite energies incudes a metallic phase in disordered AGR. On 
the other hand, the conductance exponentially decays as a function of the length like 
G oc exp(— aL). It implies Anderson localization and a is the inverse of the localization 
length, a is zero for metallic states and finite for localized states depending on the 
parameters a and S. In Fig. ^c), the scaling behavior of the conductance shows that 
there also exists metal-insulator transition induced by the disorder strength of correlated 
disorder. From finite size scaling behavior, it is easy to determine the critical a cr in term 
of S and the critical 5 cr in term of a from Fig. [3]^a-c) respectively. As shown in Fig. 
Wh), a cr is 1.8~1.9 at 5 = 0.1. The dependence of conductance on the length of long 
AGR shows a cr = 1.86, which is close to but a little larger than the value from the Dirac 
equation [31]. Moreover, the conductance of Dirac fermion under disordered potential 
with long-range correlation scales as G oc L~ v , in which n is a constant and depends on 
a and 5. The discrepancy comes from the difference of tight binding model and Dirac 
equation, however both tight binding model and Dirac equation predict metal-insulator 
transition induced by the disorder with long-range correlation. Different from neutral 
AGR, the conductance oscillates around a fixed value and monotonically decay at the 
energy shift fi is 0.2 and 1.0 respectively when the length increases from 20 to 340, as 
shown in (d). 

3.2. Dependence of the conductance on the energy shift in the presence of long-range 
correlated disorder 

Fig. HJshows how the conductance varies with the energy shift at three different strengths 
of disorder. Here energy shift in unit of the nearest-neighbor hopping amplitude t is 
linear to the gate voltage applied to AGR. The size parameter for AGR is L = 20 
and M = 40, and a = 2.5. The conductance of disordered AGR is the average of 500 
samples. For comparison, the conductance of uniform AGR is also shown. The slope of 
conductance at positive energy shift is different from the counterpart at negative energy 
shift, since the electron-hole symmetry is broken due to odd- numbered ring at the AGR- 
QW interface [35]. For weak disorder, e.g., 5 = 0.1, the conductance curve is slightly 
different from that of uniform AGR, and asymmetry in electron-hole still presents. As 
the disorder strength S increases to 0.5, the conductance decreases for most energy shift 
while increases around zero energy shift. The variance of conductance increases with the 
disorder strength. When 5 is 1.0, only the conductance near zero energy shift is highest, 
which is completely different from those above curves of the conductance versus the 
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Figure 4. Dependence of the conductance on the energy shift when 5 changes from 
to 1.0. a— 2.5. The length and width of sample is L — 20 and M — 40 respectively. 




Figure 5. Dependence of the conductance in AGR under several energy shift on the 
disorder strength S when a is 1.0, 2.0 and 2.5 respectively. The length and width of 
sample is L = 20 and M — 40 respectively. 



Fig. [5] shows the dependence of the conductance G at three different energy shifts 
on the disorder strength 5. For comparison, a =1.0, 2.0 and 2.5 are chosen. As shown 
in Fig. [5]^a), the conductance of neutral AGR slightly changes with S at a = 1.0, except 
that the variance of conductance increases with S. The conductance first increases and 
then decreases as S increaes when a is larger than 2.0. This behavior is similar to 
that in the presence of long-range disorder [11] . On the contrary, the conductance first 
remains the same at weak disorder and then decay at strong disorder in the presence 
of short-range disorder such as Anderson disorder. When the energy shift is 0.2, the 
conductance decreases as 5 increases for lowly correlated disorder, and the behavior of 
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the conductance in the presence of highly correlated disorder is similar to that of neutral 
AGR , as shown in Fig. Efb). When the energy shift is 1.0, the conductance decreases as 
the disorder strength increases for all a as shown in Fig. [5](c) . This is due to Anderson 
localization in heavily dopped AGR. 

It is well known that edge disorder in graphene induces insulating behavior in 
graphene [H]. Different from complex edge disorder, a simple edge disordered graphene 
nanoribbon is constructed, in which 5% carbon atoms (i.e. 0.05L and 0.05L) are added 
randomly to carbon atoms at the top and bottom edges respectively. Then we analyze 
the conductance of 100 such random samples in the presence of long-range correlated 
disorder with a = 2.5. For each random sample, the conductance is averaged between 
500 realizations of long-range correlation disorder, as we did for the case without edge 
disorder. On one hand, the conductance increases with the disorder strength 5(< 0.5), 
as shown in Fig. El On the other hand, the conductance decreases as the length 
increases for all 5, which no derealization occurs in the presence of highly long-range 
correlated disorder. This is contrary to the fact that a = 2.0 induces derealization 
in long-range correlated disordered graphene nanoribbon without edge disorder. Higher 
percent edge disorder results in lower conductance, while the dependence of conductance 
on 5 is similar. It indicates that no derealization occurs in edge-disordered graphene 
nanoribbon even when the disorder is correlated, since localization is induced by edge 
disorder. 




Figure 6. Dependence of the conductance on the disorder strength 6 as the length L 
varies from 20 to 60. Full and hollow symbols represent 5% and 10% edge disorder. 
a =2. 5, and the width of sample is M = 4. 



3. 3. Dependence of the conductivity on the energy shift and on the aspect ratio in the 
presence of long-range correlated disorder 

It is well known that the conductivity of uniform and neutral armchair graphene 
nanoribbon approaches ^ as the aspect ratio Mj L is far larger than 1 (121 H3] . The 
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conductivity is denned as a = Fig. [7] shows the dependence of the conductivity on 
the energy shift in the presence of correlated disorder with a = 2.5 and 5 = 1.0. Under 
such parameters as studied in Sec. Ill (A), all electronic states are localized at the 
thermodynamic limit and the conductance decreases as the length increases as shown 
in the inset of Fig. [7J therefore the conductivity approaches to zero as the aspect ratio 
approaches to zero. However, for moderate aspect ratio, the conductivity increases at 
\fi\ < 0.75 and decreases otherwise as the length increases from 20 to 100, as shown 
in Fig. [7J The transition of conductivity in term of the length indicates that the 
conductance scales as with /3 varying from less than 1 to larger than 1, because 
the localization length modulated by the energy shift is longer than and shorter than 
the length of the system L respectively. We further find that this kind of transition of 
conductivity occurs in a series of AGRs with different width. 

Fig. [8] shows the dependence of the conductivity on the aspect ratio in the presence 
of disorder with long-range correlation with a = 2.5 and 5 = 1.0. Interestingly, there 
are two different behaviors for conductivity in term of the aspect ratio. Thereinto, the 
conductivity near zero energy shift (/x=0 and 0.2) first increases and then decreases 
while the counterpart at fi=1.0 monotonically decreases, as the aspect ratio approaches 
to zero. For comparison, the corresponding behavior of the conductivity for uniform 
AGR with M = 40 is also shown. Since the conductance remains finite and constant 
for metallic AGR, the conductivity scales as 0(L/M) as M/L approaches to zero for 
all three different gate voltages. The conductivity corresponding to one quantized 
conductance is shown by the dotted black line. It is easy to find that the asymptotical 
conductivity of disordered AGR with long-range correlation is different from metallic 
uniform AGR, since localization occurs at the thermodynamic limit. On the other 
hand, the conductivity of disordered AGR at /i=1.0 decreases slower than that of 
neutral AGR with M = 39, which is semiconducting without impurity or disorder, 
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Figure 8. Dependence of the conductivity on the aspect ratio M/L. a=2.5 and 
(5=1.0. The energy shift /x is (square), 0.2 (circle) and 1.0 (upper triangle). The 
width of sample is 20 (hollow symbols) and 40 (filled symbols). Square, circle and 
upper triangle with cross stand for the conductivity of uniform AGR with M = 40. 
Half- filled lower triangle stands for uniform neutral AGR with M = 39. The dotted 
black line is the conductivity corresponding to one quantized conductance. 



while their asymptotical behavior of the conductivity should be the same because of 
localization. Compared with pure AGR, the conductivity of disordered AGR with long- 
range correlation depends smoothly on the aspect ratio in general. 

4. Summary 

Transport in disordered armchair graphene nanoribbons (AGR) with long-range 
correlation between quantum wire contact is investigated based on tight binding model. 
Metal-insulator transition is induced by disorder in neutral AGR. Thereinto, the 
conductance is one conductance quantum for metallic phase and exponentially decays 
otherwise, when the length approaches to infinity and is far longer than the width. 
Similar to the case of long-range disorder, the conductance of neutral AGR first increases 
and then decreases while the conductance of doped AGR monotonically decreases, as 
the disorder strength increases. In the presence of strong disorder, the conductivity 
depends monotonically and non-monotonically on the aspect ratio for heavily doped 
and slightly doped AGR respectively. For edge disordered graphene nanoribbon, the 
conductance increases with the disorder strength of long-range correlated disordered 
while no derealization occurs, since localization is induced by the edge disorder. 

Acknowledgment s 



This work was supported by the Research Foundation from Ministry of Education 
of China (2009EDU309002), the National Basic Research Program of China 



Effect of disorder with long-range correlation on transport in graphene nanoribbon 11 



(2012CB921704), and the National Natural Science Foundation of China under Grant 
No 11174363. 



References 

Novoselov K S et al 2004 Science 306 666 
Berger C et al 2004 J. Phys. Chem. B 108 19912 
Novoselov K S et al 2005 Nature (London) 438 197 
Zhang Y et al 2005 Nature (London) 438 201 
Seol J H, Jo I, Moore A L et al 2010 Science 328 213 

Hubert B H, Pablo J-H, Jeroen B O, Lieven M K V and Alberto F M 2007 Nature 446 56-59 
Cresti A, Nemec N, Biel B et al 2008 Nano Res. 1 361-394 

Castro Neto A H, Guinea F, Peres N M R et al 2009 Rev. Mod. Phys. 81 109-162 
Mucciolo E R, Lcwcnkopf C H 2010 J. Phys.: Condensed Matt. 22 2732201 

Abergel DSL, Apalkov V, Berashevich J, Ziegler K and Chakraborty T 2010 Adv. Phys. 59 
261-482 

Rycerz A, Tworzydlo J and Beennakker C W J 2007 Europhys. Lett. 79 57003 
Zhang Y Y, Hu J P, Bernevig B A, Wang X R, Xie X C and Liu W M 2009 Phys. Rev. Lett. 102 
106401 

Xiong S J and Xiong Y 2007 Phys. Rev. B 76 214204 
Qin Z J and Zhang G P 2011 J. Phys.: Condens. Matter 23 295301 
Hwang E H, Adam S and Das Sarma S 2007 Phys. Rev. Lett. 98 186806 

Chen J H, Jang C, Adam S, Fuhrcr M S, Williams E D and Ishigami M 2008 Nat. Phys. 4 377-381 
Martin J, Akerman N, Ulbricht G et al 2008 Nat. Phys. 4 144-148 
Anderson P W 1958 Phys. Rev. 109 1492 

Abrahams E, Anderson P W, Licciardello D C and Ramakrishnan T V 1979 Phys. Rev. Lett. 42 
673 

Wu H-L, and Philips P W 1991 Phys. Rev. Lett. 66 1366 
Philips P W and Wu H-L 1991 Science 252 1805 
Chen Xiaoshuang and Xiong S J 1993 Phys. Lett. A 179 217 
dc Moura F A B F and Lyra M L 1998 Phys. Rev. Lett. 81 3735 
Izrailev F M and Erokhin A A 1999 Phys. Rev. Lett. 82 4062 
Zhang G P and Xiong S J 2002 Euro. Phys. J. B 29 491 
Liu W S, Chen T and Xiong S J 1999 J. Phys.: Condens. Matter 11 6883 
Liu W S, Liu S Y and Lei X L 2003 Euro. Phys. J. B 33 293 
Izrailev F M, Erokhin A A and Makarov N M 2012 Physics Report 512 125-254 
Li Q Z, Hwang E H, Rossi E and Das Sarma S 2011 Phys. Rev. Lett. 107 156601 
Yan J, and Fuhrer M S 2011 Phys. Rev. Lett. 107 206601 

Cheraghchi H, Irani A H, Fazeli S M and Asgari R 2011 Phys. Rev. B 83 235430 
Federenko A A, Carpentier D and Orignac E 2012 Phys. Rev. B 85 125437 
Schomerus H 2007 Phys. Rev. B 76 045433 
Zhang G P and Qin Z J 2010 Phys. Lett. A 374 4140 
Zhang G P and Qin Z J 2011 Chem. Phys. Lett. 516 225 
Robinson J P and Schomerus H 2007 Phys. Rev. B 76 115430 
Zhang G P, Wang C Z and Ho E M 2011 Phys. Lett. A 375 1043 
Yin Y and Xiong S J 2003 Phys. Lett. A 317 507 
Son Y W, Cohen M L and Louie S G 2006 Phys. Rev. Lett. 97 216803 

Yao Y X, Wang C Z, Zhang G P, Ji M and Ho E M 2009 J. Phys.: Condens. Matter 21 235501 
Louis E, Verges J A, Guinea F, and Chiappe G 2007 Phys. Rev. B 75 085440; Evaldsson M, 
Zozoulenko I V, Xu H Y, and Heinzel T 2008 Phys. Rev. B 78 161407(R); Cresti A, and Roche 
S 2009 New J. Phys. 11 095004; Gunlycke D, and White C T 2010 Phys. Rev. B 81 075434 



Effect of disorder with long-range correlation on transport in graphene nanoribbon 12 
[42] Katsnclson M I 2006 Euro. Phys. J. B 51 157 

[43] Tworzydlo J, Trauzettcl B, Titov M, Ryccrz A and Beenakker C W J 2006 Phys. Rev. Lett. 96 
246802 



